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Abstract We theoretically investigate a Bose-condensed exciton gas out of equilibrium. 
Within the framework of the combined BCS-Leggett strong-coupling theory with the non¬ 
equilibrium Keldysh formalism, we show how the Bose-Einstein condensation (BEC) of 
excitons is suppressed to eventually disappear, when the system is in the non-equilibrium 
steady state. The supply of electrons and holes from the bath is shown to induce quasi¬ 
particle excitations, leading to the partial occupation of the upper branch of Bogoliubov 
single-particle excitation spectrum. We also discuss how this quasi-particle induction is re¬ 
lated to the suppression of exciton BEC, as well as the stability of the steady state. 

PACS numbers: 71.35.-y,03.75.Ss, 71.36.J-c 


1 Introduction 

Since the prediction of an exciton Bose-Einstein condensation (BEC) in a semiconductor— i^, 
this electron-hole pair condensate has attracted much attention as an analogous phenomenon 
to metallic superconductivity, where electron-electron Cooper-pairs play the central role. 
Although the exciton BEC has not been realized yet, recent experiments have explored sub- 
Kelvin temperature s . Thus, the realization of an exciton BEC is very promising. 

Once this Fermi condensate is realized, it is expected that one can examine various 
physical properties of this system from the weak-coupling regime to the strong-coupling 
limit, by adjusting the exciton density. At a glance, this advantage is similar to the case 
of a superfluid Fermi atomic gas, where the interaction between Fermi atoms is also tun¬ 
able by adjusting the threshold energy of a Feshbach resonance-^. In the latter, the so-called 
BCS-BEC crossover has been realized-*^, where the character of superfluidity continuously 
changes from the weak-coupling BCS type to the BEC of tightly bound molecules, with 
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increasing the strength of a pairing interactio nal***!** . However, while the cold atom system 
is usually in the equilibrium state, an exciton gas is essentially in the non-equilibrium state, 
because one always needs to continue supplying electrons and holes to the system, in order 
to compensate the decay of excitons into photonsi^. Thus, the realization of an exciton BEC 
would provide a unique opportunity to examine the BCS-BEC crossover phenomena in the 
non-equilibrium case. Since non-equilibrium properties of a Eermi condensate has recently 
been discussed in various systems, such as an exciton-polariton gas in a microcavity—, as 
well as an ultracold Fermi ga s *^* 1 *^ , an exciton BEC would also contribute to the study in 
these fields. 

In this paper, we investigate non-equilibrium properties of an exciton condensate. We 
employ a model for an electron-hole gas, that has attractive interactions between the species 
to promote pairing, and also decay to a vacuum and pumping from a bath of free fermions. 
The tunneling of particles to the vacuum effectively describes the decay of excitons in this 
model. This leakage is compensated by the supply of particles from the bath, leading to 
the steady state. In this model, we examine strong-coupling effects within the BCS-Leggett 
theory at zero bath temperature, T], = 0. Effects of the non-equilibrium steady state are also 
taken into account by using the Keldysh Green’s functio n *^i*^'*^ . In the non-equilibrium 
steady state, we examine how the exciton BEC is suppressed, to eventually disappear. As 
a signature of this suppression, we show that partial occupation of the upper branch of 
Bogoliubov single-particle excitations occurs. Throughout this paper, we take Ti = = \, 

and system volume is set to unity, for simplicity. 


2 Non-equilibrium BCS-Leggett theory in the presenee of pumping and decay 


We consider a model electron-hole gas in the BEC state, described by the Hamiltonian, 

H = H,+H,^,+Hu ( 1 ) 


where, 

H, = Y^'y;[epT^-A{t)T+-A*[t)T^]%-UY^p+pZ,, (2) 

P 9 

+ J^d>;+e;T3d>;, (3) 

p p 

Hi= E EE +h.c.j . (4) 

A=b,vP,9 i 

The electron-hole gas in the exciton-BEC phase is described by Hg in Eq. where 'Fp = 
*® ^ Nambu field, consisting the electron annihilation operator (flp.e) and the 
hole creation operator {a^_p [,). These particles are assumed to have the same mass m, as 
well as the same kinetic energy Ep = p^/{2m). T,- (i = 1,2,3) are Pauli matrices acting 
on electron-hole space, and where T± = [ti ±!T2]/2. The attractive 

interaction between an electron and a hole is modeled by a contact interaction with the 
coupling constant —U (< 0) (although the real interaction is, of course, the long-range 
Coulomb interaction). For simplicity, we ignore the repulsive interaction between electrons, 
as well as the interaction between holes. The exciton-BEC state is characterized by the order 
parameter A (t) = U'£k{a-k,h{f)‘^k.eit))- 

The model exciton-BEC gas is coupled to a bath and a vacuum described by the first 
and second term in Hem in Eq. respectively. In. Eq. ([3]l, the Nambu field = 
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Fig. 1 (Color online) Diagrammatic expressions for the self-energy, (a) (b) The solid line and 

the dashed line represent the Keldysh Green’s function in the exciton-BEC system and that in the bath or 
vacuum, respectively. The dotted line describes the electron-hole interaction —U, and solid circle represents 
the transfer matrix element Ib,v between the exciton-BEC system and the environment (consisting of a bath 
and a vacuum). 


consists of the electron annihilation operator and the hole creation op¬ 
erator in the bath (vacuum). Electrons and holes in the bath (vacuum) are assumed to 

have the same kinetic energy for simplicity. 

The momentum-independent transfer matrix element (i^) in Eq. (01 represents the 
coupling between the exciton-BEC system and the bath (vacuum). Here, we assume that 
electrons and holes of the exciton BEC at position r, tunnels to position /?,■ in the bath (vac¬ 
uum). While the leakage of particles from the exciton-BEC gas into the vacuum effectively 
describes the decay of excitons, the tunneling from the bath to the exciton-BEC works as 
the pumping, to compensate the leakage of particles into the vacuum. 

For the time dependence of the BEC order parameter A (?) in the non-equilibrium steady 
state, we employ the ansatziiiS 

A{t) = , (5) 

where Ao is taken to be real. In this case, one can formally eliminate the time depen¬ 
dence from the model Hamiltonian in Eq. (|T]| by the gauge transformation (ap,Cp,cp = 
{ap,c^,Cp)e^'^‘ . The resulting Hamiltonian has the same form as Eq. (|T]l, where Ep, Ep''', 
and A{t), are replaced by ^p = Ep — fi, ^p'' = Ep'' — )i, and 4o. respectively. 

As usual, we assume that the bath and the vacuum are huge compared to the exciton- 
BEC system, so that they are still in the thermal equilibrium state, even when they are 
coupled to the exciton-BEC system. In particular, we consider the case where the bath is 
at 7], = 0. In this case, the electron and hole distribution in the bath is simply given by the 
ordinary Fermi distribution function at Tf, = 0, /b(ffl) = S(co — [)ib “ M])’ where jUb is the 
Fermi energy in the bath. For the vacuum, since the particles are absent there, we take the 
vanishing distribution in the vacuum as /v(ffl) = 0. 

To systematically examine non-equilibrium effects on the strong-coupling exciton BEC 
at 7b = 0, it is convenient to reformulate the BCS-Leggett theory using the Keldysh Green’s 
function G„ ®)— , which obeys the Dyson equation 

Ga,j3(p,C0) = G“ „,(/>,(»)+ G° ij(p,w)Lpjji(p,CO)Gpra'(p,(o), (6) 


where the lowest-order Keldysh Green’s function G° (p, a) has the form, under the Nambu 
representation. 


G«(p.a)) = {G°}„,a'=( '■ 




CO id — 

0 


co-iS-^pTs 




(7) 
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In Eq. 0, 5 is an infinitesimally small positive number, and G^, G®, and G^, are the 
retarded, advanced, and Keldysh components, respectively. The self-energy „/(p, co) in 
Eq. ® involves effects of the electron-hole interaction (—G), as well as the coupling to the 
bath (Tb) and the vacuum (1^). In the BCS-Leggett theory, the former is treated within the 
Hartree-Fock-Bogoliubov (HFB) approximation, where the HFB self-energy E™^,{p,co) is 
diagrammatically described as Fig.[TJa), which givesiS 


E»^,(p,co) = iU'^f 

n' J 


da' 

1m 


E E 

I3,p's,s’=z 




lrlT,G(i,,j,(p',a')]T,, 


( 8 ) 


Here, = (^a,a'^p.-p' + ^a,-a'^p,p')/^’ where —a' means the opposite component 

to a'. We also include the couplings to the bath (Tb) and the vacuum (i^) in the second 
Bom approximation shown diagrammatically in Fig0b). The expression for this self-energy 
correction Z®™ is given by, after taking the random average over the tunneling position r,-. 


y^env 

^a,a 


,(p,a} 




-iyx -2iTiYA{(OT3) 
0 iyi 


■ (9) 

a,a' 


Here, {q, a) are the non-interacting Keldysh Green’s functions in the bath (X = b) and 
the vacuum {X =y), = I — 2fx{a), and 

Y)^ = nNiP^\rx\^ (10) 

describes pumping and decay effects by the coupling to the bath {X = b) and vacuum (A = v) 
(A^t is the number of tunneling positions). In obtaining Eq. ( llQb . the bath and the vacuum are 
assumed to be white, that is, each part has the energy-independent single-particle density of 
states Pi. 

Using the self-energy E = _|_£env^ obtains the self-consistent equation from 

Gii(p, a) in Eq. ([6]l, which corresponds to the BCS gap equation in the equilibrium state, 
as 

1 uV f da F(a)[a + ^p + iY\-F{-a)[a-^p-iY\ 

271 [{co-Ep)^ + Ym{a + Ep)^ + Y^ ’ ^ ^ 


where 7 = 'll, + 7v, and F{a) = %Fb{a) -I- YvFy{a). Ep = + Aq is the ordinary the 

Bogoliubov single-particle excitation spectrum. Since Eq. ill lb involves the ultraviolet di¬ 
vergence, as in the ordinary BCS gap equation, we need to renormalize the theory to elimi¬ 
nate this singularity. This is conveniently achieved by measuring the interaction strength in 
terms of the j-wave scattering length a,, which is related to the contact interaction — G as 
AKUsIm = —G/[l —UYfp{^/{2£p))]- In this scale, < 0 and 0 > ' repre¬ 

sent the weak-coupling side and the strong-coupling side, respectively (where kp the Fermi 
momentum), pc is a cutoff momentum. 

Following the BCS-Leggett theory in the equilibrium stat c^i**^ , we solve the “gap equa¬ 
tion” l[II]l, together with the equation for the total number N of electrons and holes, which 
is obtained from (1,1)-component (in Nambu space) of the lesser Green’s function G^ = 
—Gi 1 -|- G 22 + Gi 2 as— 


N = -2i 




n 


da [{a + ^pf + Al + Y^\yhfh{oi)+^l[yh{'^-fh{-oi)) + 7^ 
K [{a-EpY + Y^][{a + EpY + Yl 


7Z 


(12) 
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Fig. 2 (Color online) (a) Exciton-BEC order parameter Aq in the unitarity limit ((kpas) * = 0) at Tb = 0, as 
a function of the decay parameter We take yb/ep = 10“^, where ep is the Fermi energy of a free electron- 
hole gas. The inset shows /ip and (t*) Intensity of the occupied spectral weight L(p, (o), normalized by EjT*. 
We set 7y/yj, = 0.3. The other parameters are the same as those used in panel (a). 


to self-consistently determine 2io,MB and fi, for a given parameter set (N,as,yb,'Y^). In this 
regard, we briefly note that the gap equation dl 11 1 is a complex equation, so that the real 
and imaginary components give two independent equations. Thus, together with the number 
equation {H, one may safely determine these three quantities in a consistent manner. 


3 Depairing effect in the non-equilibrium exciton-BEC phase 

Figure 12 a) shows the self-consistent solutions for the coupled equations lllU and Ill2li . In 
the equilibrium case (which is realized when = 0 ), the inset shows the expected relation 
jLib = M- When this chemical equilibrium condition is satisfled, the imaginary part of the gap 
equation O identically vanishes (Note that Fb(—®) = —Fb(®) when jtib = M-)! and the 
real part of this equation is reduced to the ordinary BCS gap equation in the limit, yt, —y 0. 

When the leakage of electrons and holes into the vacuum is turned on {y^ > 0), the bath 
must supply particles to the exciton-BEC system so as to compensate this leakage. Because 
of this, the chemical potential jUb in the bath rapidly increases when y^ is small but Unite. 
(See the region around /v = 0 in the inset in Fig. [2a).) We also see in Fig. |2a) that the 
magnitude of the exciton-BEC order parameter Aq decreases, which physically means that 
the BEC phase is suppressed in the non-equilibrium state, even when the particle loss is 
compensated by the bath. In the unitarity limit shown in Fig. |2a), the exciton-BEC phase 
disappears when 7 v/')l, > 0.96. 

The decrease of the order parameter Aq in the presence of pumping ( 7 , > 0) and decay 
(7v > 0) implies that the depairing effect occurs in the non-equilibrium exciton-BEC. To 
confirm this in a simple manner, we consider the occupied spectral weight L(p, co), which 
directly gives information about occupied single-particle states, given by. 


L(p,co) 


-iGu {P, ®) 

^/b(®)C(p,co) + 


^[l-/b(-®)] + ^|[l-C(p,®)] 


A{p,co), 


(13) 
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Here, ^(p, co) = [{co + + y^]/[(co + + ^o+ fl and 


(ffl + £'p)2 + y2 

is the spectral weight. In the equilibrium state (/v = 0), the occupied spectral weight L{p, co) 
is reduced to (Note that pb = M-) 


A(p,co) = 


2n 


1 + - 


+ ■ 


{(o-EpY + y^ In 


1 -- 


L{p,a) = &{-w)A{p,a), (15) 

where @{x) is the step function. Equation ( 115b indicates that the upper Bogoliubov single¬ 
particle band (o = Ep) is almost unoccupied, when %/Aq <C 1. However, in the non¬ 
equilibrium state. Fig. [2jb) shows the sizable occupation of this branch, when the “biased 
voltage” Pb “ M exceeds the energy gap 4o- Since the total number N of electrons and holes 
remains unchanged in the steady state, the appearance of the Bogoliubov (unpaired) quasi¬ 
particles is accompanied by the decrease of the condensate fraction in the exciton-BEC. 

The condition — fj. > Aq is also obtained as the condition to maintain the steady state, 
when %/Aq <C 1. In this limiting case, the imaginary part of Eq. (Illb gives 

= (16) 

p p 

To satisfy this, the condition 


^l^,-^^>Mm[Ep]=Ao, (17) 

is necessary, otherwise the left hand side of Eq. ( 116b vanishes identically. This condition 
is similar to the existence of threshold voltage in a normal-metal-superconductor junction, 
where quasi-particle current can flow through the junction when the biased voltage (multi¬ 
plied by the electric charge) exceeds the BCS energy gap^. When Eq. ([mi is satisfied, the 
“tunneling current” from the bath to the exciton BEC gas occurs in the energy region, 

A)<ffl<Pb-M- ( 18 ) 

As mentioned previously, to maintain the non-equilibrium steady state, the leakage of 
particles into the vacuum must be compensated by the supply of electrons and holes from 
the bath. In this regard, we note that p.}, — p. is not sensitive to the decay parameter except 
for 7v/ll, ^ 1- (See the inset in Fig.|2[a).) Thus, when y^ increases under a fixed value of 
the pumping parameter 7 ,, the exciton-BEC order parameter Aq must decrease, in order to 
increase the tunneling flow from the bath by widening the region given in Eq. (EUl- When 
7 v is very large, this compensation mechanism no longer works, leading to the vanishing 
exciton-BEC, as seen in Fig.|2[a). 


4 Summary 

To summarize, we have investigated a model exciton-BEC gas which is coupled to a bath 
and a vacuum. Reformulating the BCS-Leggett strong-coupling theory by using the Keldysh 
Green’s function, we have examined how the leakage of electrons and holes into the vacuum, 
as well as the supply of these particles from the bath, affect the stability of the exciton- 
BEC. In the non-equilibrium steady state, we showed that the BEC order parameter Aq 
decreases with increasing the leakage of electrons and holes into the vacuum, to eventually 
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vanish. We also examined the occupied spectral weight function. In the non-equilibrium 
BEC state, we found that the partial occupation of the upper branch of Bogoliubov single¬ 
particle excitations occurs, indicating the depairing of electron-hole pairs. We also pointed 
out that this phenomenon is deeply related to, not only the decrease of the exciton-BEC 
order parameter, but also the condition for the stability of the non-equilibrium BEC state. 

In this paper, we dealt with a model electron-hole mixture with a contact interaction at 
TJ, = 0. Extension of this simple treatment to include the realistic long range Coulomb inter¬ 
action, as well as pairing fluctuations at finite temperatures, is an exciting future problem. 
Since Fermi condensates out of equilibrium has recently discussed in various systems, not 
only in an electron-hole gas, but also in an exciton-polariton condensate in a microcavity, 
as well as in an ultracold Fermi gas, our results would be useful for the development of 
non-equilibrium Fermi condensates. 
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